Electron-acoustic plasma waves: 
oblique modulation and envelope solitons * 
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Theoretical and numerical studies are presented of the amplitude modulation of electron-acoustic 
waves (EAWs) propagating in space plasmas whose constituents are inertial cold electrons, Boltz- 
mann distributed hot electrons and stationary ions. Perturbations oblique to the carrier EAW prop- 
agation direction have been considered. The stability analysis, based on a nonlinear Schrodinger 
equation (NLSE), reveals that the EAW may become unstable; the stability criteria depend on 
the angle 6 between the modulation and propagation directions. Different types of localized EA 
excitations are shown to exist. 
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I. INTRODUCTION 

Electron acoustic waves (EAWs) are high-frequency (in 
comparison with the ion plasma frequency) electrostatic 
modes in plasmas where a 'minority' of inertial cold 
electrons oscillate against a dominant thermalized back- 
ground of inertialess hot electrons providing the neces- 
sary restoring force. The phase speed Vph of the EAW 
is much larger than the thermal speeds of cold elec- 
trons (vth,c) and ions (vth.i), but is much smaller than 
the thermal speed of the hot electron component {vth,h)] 
Vth,a = iTa/niaY/'^ , where a = c,h,i (ttIq denotes the 
mass of the component a; the Boltzmann constant ks is 
understood to precede the temperature Tq, everywhere). 
Thus, the ions may be regarded as a constant positive 
charge density background, providing charge neutrality. 
The EAW frequency is typically well below the cold elec- 
tron plasma frequency, since the wavelength is larger 
than the Debye length A^, = (T/i/47rn/ioe^)^^^ involving 
hot electrons (tIq. denotes the particle density of the com- 
ponent a everywhere). 

The linear properties of the EA waves are well under- 
stood [2-5]. Of particular importance is the fact that 
the EAW propagation is only possible within a restricted 
range of the parameter values, since both long and short 
wavelength EAWs are subject to strong Landau damp- 
ing due to resonance with either the hot or the cold (re- 
spectively) electron component. In general, the EAW 
group velocity scales as Vph — Vth,h \/ nc/uh] therefore, 
the condition Vth,c ^ Vph ^ Vth,h immediately leads to 
a stability criterion in the form: Tc/Th ^ \J ric/nh ^ 1. 
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A more rigorous investigation [4] reveals that EAWs will 
be heavily damped, unless the following (approximate) 
conditions are satisfied: Th/Tc ^ 10 and 0.2 ^ nc/rif. < 
0.8 (where rie = ric + nt)- Even then, however, only 
wavenumbers k between, roughly, 0.2k£)^c and 0.6 ku ^ 
(for Th/Tc — 100; see in Ref. 4b), will remain weakly 
damped [kD,c = {imicfie'^ /TcY^^ = X'^^ obviously de- 
notes the cold electron Debye wavenumber]. The sta- 
ble wavenumber value range is in principle somewhat ex- 
tended with growing temperature ratio Th/Tc, see the 
exhaustive discussion in Refs. 4 and 5. 

As far as the nonlinear aspects of EAW are concerned, 
the formation of coherent EA structures has been con- 
sidered in a one-dimensional model involving cold 0] 
or finite temperature Q ions. Furthermore, such non- 
envelope solitary structures, associated with a localized 
compression of the cold electron density, have been shown 
to exist in a magnetized plasma [9-11]. It is worth 
noting that such studies are recently encouraged by the 
observation of moving EAW-related structures, reported 
by spacecraft missions e.g. the FAST at the auroral re- 
gion [12 - 14], as weU as the GEOTAIL and POLAR 
earlier missions in the magnetosphere [15 - 17]. How- 
ever, although most solitary wave structures observed are 
positive potential waves (consistent with an electron hole 
image) , there have also been some negative potential and 
low velocity structure observations, suggesting that some 
other type of solitary waves may be present in the mag- 
netosphere [17, 18]. These structures are now believed 
to be related to EA envelope solitary waves, for instance 
due to trapping and modulation by ion acoustic density 
perturbations p. ITsI]. 

Amplitude modulation is a long-known generic feature 
of nonlinear wave propagation, resulting in higher har- 
monic generation due to nonlinear self-interactions of the 
carrier wave in the background medium. The standard 
method for studying this mechanism is a multiple space 
and time scale technique [l9l . I20I ] , which leads to a non- 
linear Schrodinger equation (NLSE) describing the evo- 
lution of the wave envelope. Under certain conditions, 
it has been shown that waves may develop a Benjamin- 
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Feir-type (modulational) instability (MI), i.e. their mod- 
ulated envelope may be unstable to external perturba- 
tions. Furthermore, the NLSE-based analysis, encoun- 
tered in a variety of physical systems [21 - 23], reveals 
the possibility of the existence of localized structures (en- 
velope solitary waves) due to the balance between the 
wave dispersion and nonlinearities. This approach has 
long been considered with respect to electrostatic plasma 
waves [20, 24 - 30]. 

In this paper, we study the occurrence of modulational 
instability as well as the existence of envelope solitary 
structures involving EAWs that propagate in an unmag- 
netized plasma composed of three distinct particle popu- 
lations: a population of 'cold' inertial electrons (mass me, 
charge — e), surrounded by an environment of 'hot' (ther- 
malized Boltzmann) electrons, moving against a fixed 
background of ions (mass m^, charge Qi = +Zie), which 
provide charge neutrality. These three plasma species 
will henceforth be denoted by c, h and i, respectively. 
By employing the reductive perturbation method and ac- 
counting for harmonic generation nonlinearities, we de- 
rive a cubic Schrodingcr equation for the modulated EA 
wave packet. It is found that the EAWs are unstable 
against oblique modulations. Conditions under which 
the modulational instability occurs are given. Possible 
stationary solutions of the nonlinear Schrodinger equa- 
tion are also presented. 



Re-scaling all variables and developing around <E> = 0, 
Eqs. O - ® can be cast in the reduced form 

9n „ , 

--KV.(nu) . 0, 
du 

— -fu-Vu = V0, 

and 

- <^ + i</)2^i03+ /3(n-l), (5) 

where all quantities are non-dimensional: n = nc/ucfi, 
u — Uc/vo and </> = $/$o; the scaling quantities are, 
respectively: the equilibrium density rico, the 'elec- 
tron acoustic speed' vq = Cs^h = {kBTh/meY^^ and 
^0 — {kBTh/e). Space and time are scaled over the De- 
bye length XD,h = (fcsT'ft/47rn;i,oe^)^/^ and the inverse 
plasma frequency uj~l = Xdji/cs = (47rn/i,oe^/TOe)~^/^, 
respectively. The dimensionless parameter (3 denotes 
the ratio of the cold to the hot electron component i.e. 
[3 = nc/rih. Recall that Landau damping in principle 
prevails on both high and low values of beta (cf. the 
discussion in the introduction). According to the results 
of Ref. [4b], for undamped EA wave propagation one 
should consider: 0.25 ^ /? ^ 4. 



II. THE MODEL EQUATIONS 



III. PERTURBATIVE ANALYSIS 



Let us consider the hydrodynamic-Poisson system of 
equations for the EAWs in an unmagnetized plasma. The 
number density Uc of cold electrons is governed by the 
continuity equation 



— + V • u, = , 
at 

where the mean velocity Uc obeys 

duc e 

— + Uc • Vuc = — V $ . 

at m„ 



(1) 



(2) 



Here, the wave potential <& is obtained from Poisson's 
equation 

= -4n '^qsUs = 4:TTe{nc + Uh ~ ZiUi) . (3) 

We assume immobile ions {rii = iiifi = constant) and a 
Boltzmann distribution for the hot electrons, i.e. w 
nil ft exp{e^/kBTii) [Th is electron temperature, ks is 
the Boltzmann constant), since the EAW frequency is 
much higher than the ion plasma frequency, and the EA 
wave phase velocity is much lower than the electron ther- 
mal speed [Th/meY^'^ . The overall quasi-neutrality con- 
dition reads 



Let S be the state (column) vector (n, u, (/))^, de- 
scribing the system's state at a given position r and 
instant t. Small deviations will be considered from 
the equihbrium state S'(°) = (1, 0, 0)^ by taking S = 

^(0) + g^d) + g2^(2) ^ ^ giO) ^ J2n=l^"S("\ 

where e ^ 1 is a smallness parameter. Following the 
standard multiple scale (reductive perturbation) tech- 
nique 0], we shall consider a set of stretched (slow) 
space and time variables C. = e{x — Xt) and r = e^t,, 
where A is to be later determined by compatibility re- 
quirements. All perturbed states depend on the fast 
scales via the phase 9i — k ■ r — ujt only, while the slow 

(n) 

scales only enter the Z— th harmonic amplitude Sf , viz. 
g{n) ^ J^Z^oo 'S'f"'(C, t) the reality condi- 

tion S^^\'' = S';''"'* is met by all state variables. Two 
directions are therefore of importance in this (three- 
dimensional) problem: the (arbitrary) propagation di- 
rection and the oblique modulation direction, defining, 
say, the a;— axis, characterized by a pitch angle 9. The 
wave vector k is thus taken to be k = (k^, ky) = 
{k cos 9, k sin6'). 

Substituting the above expressions into the system of 
equations |SJl and isolating distinct orders in e, we obtain 
the nth-order reduced equations 



nc,o + nh.o - Zi ni,Q = . 



(4) 
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(n-l) 



du 



(n-2) 



OO oo 

E E 

n' — 1 I' — — ex. 



(n-n'-l) 9u|r 



5C 



= 0, (7) 



and 



An) 



+2ilk: 



^ oo oo 

n' —1 l' — — oQ 

oo oo 

n',n" — 1 l'.l" — — QO 



4""^ - 0. (8) 



For convenience, one may consider instead of the vec- 
torial relation Q the scalar one obtained by taking its 
scalar product with the wavenumber k. 

The standard perturbation procedure now consists in 
solving in successive orders ~ e" and substituting in sub- 
sequent orders. For instance, the equations for n — 2, 
I = 1 



ilojn'p + ilk ■ up' = , 



-iluju 



(1) 



ilk(Pl 



(1) 



(9) 
(10) 



and 



(Z^/fc^ -f 1) - (3n 



,(1) 







provide the familiar EAW dispersion relation 
2 /3fc' 



fc2 + 1 



i.e. restoring dimensions 



fc2 + kl 



(13) 



where c = Cs^c/^d,c = (47rnc,oe^/we)^/^ (associated 
with the cold component), and determine the first har- 



monics of the perturbation viz. 



,(1) 



1 



k 



cos f n 



(1) 



,(1) 



u; „ = — smwnp' . (14) 
k 



Proceeding in the same manner, we obtain the second 
order quantities, namely the amplitudes of the second 
harmonics 5*2 and constant ('direct current') terms Sq , 

(2) 

as well as a non-vanishing contribution to the first 
harmonics. The lengthy expressions for these quantities, 
omitted here for brevity, are conveniently expressed in 



terms of the first-order potential correction 



a(i) 



The 



equations for n ^ 2, I ~ 1 then provide the compatibility 
condition: A = Vg{k) — = uj'{k)cos9 = ^p- cos9; A 
is, therefore, the group velocity in the x direction. 



IV. DERIVATION OF THE NONLINEAR 
SCHRODINGER EQUATION 

Proceeding to the third order in e (n — 3), the equa- 
tions for I = 1 yield an explicit compatibility condition 
in the form of the Nonlinear Schrodinger Equation 



dr 



P 



(15) 



where ip denotes the electric potential correction 
The 'slow' variables {C, t} were defined above. 

The group dispersion coefficient P is related to the 
curvature of the dispersion curve as P = ^ = 



uj"{k) cos^d +uj'{k) 



the exact form of P reads 



^ ' P2w\k 



1 



1 + 3 - cos^ d 



(16) 



The nonlinearity coefficient Q is due to the carrier wave 
self-interaction in the background plasma. Distinguish- 
ing different contributions, Q can be split into three dis- 
tinct parts, viz. Q = Qq + Qi + Q2, where 



(11) Qo 



(12) 



1 



2 /33 fc2 (1 + /j2)3 _ C0S2 e 

(1+ fc2)4(i + 2/3+ k^) 
13^ + 4/3(1 + fc2)3 + 4(1 + k^C^ + 2fc2 + k^) 

(17) 
(18) 



Q2 = - 



4/3fc2 

fc2 



12^3 



/3fc6 fc2 



+ 7^ +3(3 + 8fc2) 



(19) 
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We observe that only the first contribution Qq, related 
to self- interaction due to the zeroth harmonics, is angle- 
dependent, while the latter two - respectively due to the 
cubic and quadratic terms in - are isotropic. Also, 
Q2 is negative, while Qq, Qi are positive for all values 
of k and (3. For parallel modulation, i.e. 9 — 0, the 
simplified expressions for P\g^Q and Q\g^Q are readily 
obtained from the above formulae; note that P\g^Q < 0, 
while Qlg^Q, even though positive for fc — > (see below), 
changes sign at some critical value of k. 

A preliminary result regarding the behaviour of the co- 
efficients P and Q for long wavelengths may be obtained 
by considering the limit of small fc <C 1 in the above 
formulae. The parallel {9 — 0) and oblique {9 ^ 0) mod- 
ulation cases have to be distinguished straightaway. For 
small values of fc (fc ^ 1), P is negative and varies as 



P 



8=0 



(20) 



in the parallel modulation case (i.e. 9 = 0), thus tending 
to zero for vanishing fc, while for 9 ^ Q, P in positive and 
goes to infinity as 



9^0 



VP . 2, 

TTT sm ( 
2fc 



(21) 



for vanishing fc. Therefore, the slightest deviation by 9 
of the amplitude variation direction with respect to the 
wave propagation direction results in a change in sign of 
the group-velocity dispersion coefficient P. On the other 
hand, Q varies as ~ 1/fc for small fc ^ 1. For 9 0, Q 
is negative 



Q\ 



¥0 



(22) 



while for vanishing 9, the approximate expression for Q 
changes sign, i.e. 



1 9=0 



(23) 



In conclusion, both P and Q change sign when 'switch- 
ing on' theta. Since the wave's (linear) stability profile, 
expected to be infiuence by obliqueness in modulation, 
essentially relies on (the sign of) the product PQ (see 
below), we see that long wavelengths will always be sta- 
ble. 



V. STABILITY ANALYSIS 

The standard stability analysis [20, 21, 31] consists in 
linearizing around the monochromatic (Stokes's wave) so- 
lution of the NLSE 1/; = 7/. e^'^l'^l'"' + c.c. , (notice 
the amplitude dependence of the frequency) by setting 
ip = ■00 + e "01 J and taking the perturbation 0i to be 
of the form: 01 = -01,0 e**^''''^"'^) -|- c.c. (the perturba- 
tion wavenumber fc and the frequency lo should be dis- 
tinguished from their carrier wave homologue quantities. 



denoted by fc and oj). Substituting into (|15|l . one thus 
readily obtains the nonlinear dispersion relation 



= p'k' e - 



2|l0oP 



(24) 



The wave will obviously be stable if the product PQ 
is negative. However, for positive PQ > 0, insta- 
bility sets in for wavenumbers below a critical value 

kcr — ^2p|-0o|, i.e. for wavelengths above a thresh- 
old: Xcr = 27r/fccr; defining the instability growth rate 
(T = \Imuj{k)\, we see that it reaches its maximum value 
for fc = kcr/V^, viz. 



(Jmax = \lmuj\^^ 



72 = IQII0OI 



(25) 



We conclude that the instability condition depends only 
on the sign of the product PQ, which may now be stud- 
ied numerically, relying on the exact expressions derived 
above. 

In figures d to El we have depicted the PQ = bound- 
ary curve against the normalized wavenumber k/kz) (in 
abscissa) and angle 9 (between and tt) ; the area in black 
(white) represents the region in the (fc — 9) plane where 
the product is negative (positive), i.e. where the wave 
is stable (unstable). For illustration purposes, we have 
considered w a wide range of values of the wavenumber 
fc (normalized by the Debye wavenumber k^ji, neverthe- 
less, recall that the analysis is rigorously valid in a quite 
restricted region of (low) values of fc. Modulation angle 9 
is allowed to vary between zero and tt (see that all plots 
are ^- periodic). 

As analytically predicted above, the product PQ is 
negative for small fc, for all values of theta; long wave- 
lengths will always be stable. The product possesses pos- 
itive values for angle values between zero and 9 ~ 1 rad 
~ 57°; instability sets in above a wavenumber thresh- 
old which, even though unrealistically high for 9 = 0, is 
clearly seen to decrease as the modulation pitch angle 
9 increases from zero to approximately 30 degrees, and 
then increases again up to « 57°. Nevertheless, beyond 
that value (and up to 7r/2) the wave remains stable; this 
is even true for the wavenumber regions where the wave 
would be unstable to a parallel modulation. The inverse 
effect is also present: even though certain fc values corre- 
spond to stability for 9 = Q, the same modes may become 
unstable when subject to an oblique modulation (0^0). 
In all cases, the wave appears to be globally stable to 
large angle 9 modulation (between 1 and 7r/2 radians, 
i.e. 57° to 90°). 

It is interesting to trace the influence of the percentage 
of the cold electron population (related io (3 = ric/nh) on 
the qualitative remarks of the preceding paragraph. For 
values of (3 below unity, there seems to be only a small 
effect on the wave's stability, as described above; cf. figs, 
n 121 As a matter of fact, (3 < \ appears to be valid in 
most reports of satellite observations, carried out at the 
edges of the Auroral Kilometric Radiation (AKR) region 
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(where the hot and cold electron population co-existence 
is observed) 0, 0, Q ; furthermore, theoretical studies 
have suggested that a low fi value (in addition to a high 
hot to cold electron temperature ratio) are conditions 
ensuring EAW stability i.e. resistance to damping 
Nevertheless, notice for rigor that allowing for a high frac- 
tion of cold electrons (/? ^ 3.5) leads to a strong modi- 
fication of the the EA wave's stability profile, and even 
produces instability in otherwise stable regions; cf. fig. 
13 (where an unrealistic value of /9 = 5 was considered). 
In a qualitative manner, adding cold electrons seems to 
favour stability to quasi-parallel modulation (small ^), 
yet allows for instability to higher 9 oblique modulation; 
cf. figs. ntoOl Since the black/white regions in the fig- 
ures correspond to dark/bright type solitons (see below), 
we qualitatively deduce that a solitary wave of either type 
may become unstable in case of an important increase in 
minority electron component, i.e. well above /3 = 1; see 
fig. El Notice that the critical value of the cold-to-hot 
electron number ratio (3 in order for such phenomena to 
occur may be quite low if oblique modulation is consid- 
ered; see e.g. fig. Eb. 



VI. ENVELOPE SOLITARY WAVES 

The NLSE H15() is known to possess distinct types of 
localized constant profile (solitary wave) solutions, de- 
pending on the sign of the product PQ. Following Rcf. 
|3lL l3^ . we seek a solution of Eq. H15|) in the form 
iP{(,t) = a//9(C,t) e'^'^'"^), where p, a are real variables 
which are determined by substituting into the NLSE and 
separating real and imaginary parts. The different types 
of solution thus obtained are summarized in the follow- 
ing- 

For PQ > we find the (bright) envelope soliton [33 

p = p, sech' (^^) , - ^ CH^+lu')r] , 

representing a bell - shaped localized pulse travelling at 
a speed u and oscillating at a frequency f2 (for u = 0). 
The pulse width L depends on the (constant) maximum 
amplitude square po as 



L = 



2P 



'Po 



(27) 



For PQ < we have the dark envelope soliton (hole) 



P = Pi 



^-2P 



1 — sech 



uC 



21 C-UT 



U 



= pi tank 



21 C-UT 



- 2PQpi 1 T 



(28) 



representing a localized region of negative wave density 
(shock) travelling at a speed u; this cavity traps the 



electron-wave envelope, whose intensity is now rarefac- 
tive i.e. a propagating hole in the center and constant 
elsewhere. Again, the pulse width depends on the maxi- 
mum amplitude square pi via 



L' 



P 



Qpi 



(29) 



Finally, looking for velocity-dependent amplitude solu- 
tions, for PQ < 0, one obtains the grey envelope solitary 
wave ^ 



P 

e 



P2 
1 

2P 



1 — sech^ I — 



-S sin-i 



- 2PQp2j T + eio 
a tanh( ''^"^ ) 
secK^ { ^-jjr- 



(30) 



which also represents a localized region of negative wave 
density; Oio is a constant phase; S denotes the product 
S = sign P x sign [u ~ Vq). In comparison to the dark 
soliton (|28|l . note that apart from the maximum ampli- 
tude p2, which is now finite (i.e. non-zero) everywhere, 
the pulse width of this grey-type excitation 



P 



QP2 

now also depends on a, given by 



1 



1 



2PQ p2 



in' - Vi) < 1 



(31) 



(32) 



(PQ < 0), an independent parameter representing the 
modulation depth (0 < a < 1). Vq is an independent real 
constant which satisfies the condition 

Vo - V^\PQ\P2 <u<Vo + V2\PQ\p2 ; 

for Vq = u, we have a = 1 and thus recover the dark 
soliton presented in the previous paragraph. 

Summarizing, we see that the regions depicted in figs, 
actually also distinguish the regions where different 
types of localized solutions may exist: bright (dark or 
grey) solitons will occur in white (black) regions (the dif- 
ferent types of NLS excitations are exhaustively reviewed 
in 32]). Soliton characteristics will depend on dispersion 
and nonlinearity via the P and Q coefficients; in particu- 
lar, the sign/absolute value of the ratio P/Q provides, 
as we saw, the type (bright or dark-grey)/width, re- 
spectively, of the localized excitation. Therefore, regions 
with higher values of \P/Q\ will support wider (spatially 
more extended) localized excitations of either bright or 
dark/grey type - see fig. 0] Solitons of the latter type 
(holes) appear to be predominant in the long wavelength 
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region which is of interest here, in agreement with obser- 
vations, yet may become unstable and give their place to 
(bright) pulses, in the presence of oblique perturbation 
(figs, no and/or local ndrih- value irregularities (fig. 
O. In the short wavelength region, these qualitative re- 
sults may still be valid, yet quantitatively appear to be 
rather questionable, since the wave stability cannot be 
taken for granted due to electron Landau damping. Nev- 
ertheless, even so, the EAWs are known to be less heavily 
damped than Langmuir waves , and may dominate the 
space plasma (high) frequency spectrum in the presence 
of different temperature electron populations. 

VII. CONCLUSIONS 

This work has been devoted to the study of the modu- 
lation of EAWs propagating in an unmagnetized space 
plasma. Allowing for the modulation to occur in an 
oblique manner, we have shown that the conditions for 
the modulational instability depend on the angle between 
the EAW propagation and modulation directions. In 
fact, the region of parameter values where instability oc- 
curs is rather extended for angle 9 values up to a certain 
threshold, and, on the contrary, smeared out for higher 
9 values (and up to 90 degrees, then going on in a ^ - 
periodic fashion). 

Furthermore, we have studied the possibility for the 
formation of localized structures (envelope EAW solitary 
waves) in our two electron system. Distinct types of lo- 
calized excitations (envelope solitons) have been shown 
to exist. Their type and propagation characteristics de- 
pend on the carrier wave wavenumber k and the mod- 



ulation angle 9. The dominant localized mode at long 
wavelengths appears to be a rarefactive region of negative 
wave intensity (hole), which may however become unsta- 
ble to oblique modulation or variations of the nd ratio. 
It should be mentioned that both bright and dark/grey 
envelope excitations are possible within this model; thus, 
even though the latter appear to be rather privileged 
within the parameter range where waves are expected 
not to be heavily damped, the former may exist due 
to oblique amplitude perturbations. In conclusion, we 
stress that the qualitative aspects of the observed en- 
velope solitary structures are recovered from our simple 
fluid model. The present investigation can be readily ex- 
tended to include the effects of the geomagnetic field, a 
tenuous electron beam, and on dynamics on the ampli- 
tude modulation of the EAWs. The magnetic field effects 
would produce three-dimensional NLSE in which the lon- 
gitudinal and transverse (to the external magnetic field 
direction) group dispersions would be different due to the 
cold electron polarization effect. The harmonic genera- 
tion nonlinearities would also be modified by the presence 
of the external magnetic field. 
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Figure Captions 

Figure 1. 

The product PQ = contour is depicted against the 
normalized wavenumber k/ku (in abscissa) and angle 9 
(between and tt); black (white) colour represents the 
region where the product is negative (positive), i.e. the 
region of linear stability (instability). Furthermore, black 
(white) regions may support dark (bright)-type solitary 
excitations. This plot refers to a realistic cold to hot 
electron ratio equal to (3 = 0.5 (i.e. one third of the 
electrons are cold). 

Figure 2. 

Similar to fig. ^ for (3 — 1. 
Figure 3. 

Similar to figures ^ |2] considering a very strong pres- 
ence of cold electrons {j3 = 5). Notice the appearance 
of instability (bright) regions for large angle values and 
long wavelengths. 

Figure 4. 

Contours of the ratio P/Q - whose absolute value is 
related to the square of the soliton width, see (|77|) , 
- are represented against the normalized wavenumber 
k/ko.h and angle 9. See that the negative values cor- 
respond to two branches (lower half), so that the varia- 
tion of P/Q, for a given wavenumber k, does not depend 
monotonically on 6*. P — 0.5 in this plot. 

Figure 5. 

The P/Q coefficient ratio, whose sign/absolute value 
is related to the type/width of solitary excitations, is de- 
picted against the cold-to-hot electron density ratio (3. 
The wavenumber is chosen as k/ku.h = 0.7. (a) 6* = 0° 
(parallel modulation): only dark- type excitations exist 
{PQ < 0); their width increases with /3. (b) 9 = 60° 
(oblique modulation): bright/dark excitations exist for 
(3 below/above /3cr ~ 0.8. The bright/dark soliton width 
increases/decreases with /?. (c) 9 = 90° (transverse mod- 
ulation). A rather (unacceptably) high value of f3 was 
taken, to stress the omnipresence of dark-type excita- 
tions. 
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